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OBJECT:  The  purpose  of  this  report  is  to  present  an  investigation 
•f  one  of  the  classic  papers  on  one-dimensional  shock  waves 
in  order  to  promote  understanding  of  this  phenomenon,  and  to 
present  also  an  investigation  of  one  of  the  modern  theories 
of  propagation  of  spherical  shock  waves. 

SUMMARY:  R.  Becker's  article  "Stosswelle  and  Detonation"  (1921)  is 
examined  in  detail  with  respect  to  the  theory  of  shock 
waves.  The  Hugoniot  relations  are  developed  and  are 
shown  to  be  independent  of  friction  and  heat  transfer 
effects  in  the  shock  front. 

The  theory  of  propagation  of  spherical  shocks,  as  deve¬ 
loped  by  Kirkwood  and  Brinkley  (1947)  is  discussed. 
Comparison  of  theoretical  and  experimentally  measured 
quantities  is  shown  to  demonstrate  the  validity  of  the 
theory. 

A  glossary  of  hydrodynamical  equations  and  a  detailed 
annotated  bibliography  is  also  given. 
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Section  I 


ONE  DIMENSIONAL  SHOCK  WAVES 

INTRODUCTION 

It  seems  that  Poisson1  was  the  first  to  find  a  simple  wave  solution  of 
the  differential  equation  of  one-dimensional  flow  of  a  gas  at  constant 
temperature,  and  J.  Challis^  observed  that  such  a  solution  equation  can 
not  always  be  solved  uniquely  for  the  velocity  u. 

To  remedy  this  situation,  Stokes^  proposed  to  assume  a  discontinuity  in 
the  velocity.  He  used  the  laws  of  conservation  of  mass  and  momentum  to 
deduce  two  discontinuity  conditions  for  an  isothermal  gas.  Stokes  wa3 
the  first  to  suggest  the  possibility  that  discontinuous  pressure  waves 
might  be  propagated  with  velocities  greater  than  that  of  sound.  He 
argued  further  that  discontinuities  would  never  occur  physically,  be¬ 
cause  any  tendency  to  form  such  a  discontinuity  would  be  counteracted 
by  viscous  forces, 

S.  Earnshaw4  first  developed  the  laws  of  propagation  of  waves  of  finite 
amplitude  (as  distinguished  from  sound  waves,  which  are  assumed  to  be 
infinites  sizcal)  and  investigated  irath  erratic  ally  the  building  up  of  a 
discontinuity.  He  realized  from  investigation  of  the  flow  of  gases  for 
which  pressure  i3  a  function  of  density  only,  that  since  the  local  velo¬ 
city  of  propagation  increases  across  a  compression  wave,  such  a  wave 
would  always  be  "gaining”  on  its  front  until  a  discontinuity  would  form, 

1  S*  D.  Poisson  -  Memoire  sur  la  theorie  du  son.  Journal  de  l'ecole 
polytechnique,  14^me  cahier,  7,  319-392  (1808) 

2  J.  Challis  -  On  the  velocity  of  sound.  Phil-Mag,  32,  494-499  (1848) 

3  G.  G.  Stokes  -  On  a  difficulty  in  the  theory  of  sound.  Phil-Mag. 

33,  349-356  (1848) 

4  S.  Earnshaw  -  On  the  Mathematical  Theory  of  Sound.  Phil-Trans. 

A  150,  133-148  (1850) 
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B.  RiemannS,  without  knowledge  of  Eamsh&w's  work,  developed  his  own 
theory  of  the  simple  wave,  and  obtained  the  general  solution  to  the 
flow  problem  by  introducing  what  is  now  known  as  "Riemann  invariants.” 

He  rediscovered  the  theory  of  shocks  and  elaborated  on  it,  but_i 
rectly  assumed  the  transition  across  a  shock  to  be  isentropic. 

W.  J.  M.  Rankine^*  developed  conditions  for  conservation  of  mass,  momen¬ 
tum,  and  energy  across  the  discontinuity.  He  showed  that  no  steady 
adiabatic  process  in  which  the  only  forces  are  pressure-forces  can 
represent  a  discontinuous  change  over  a  small  finite  region  from  one 
constant  state  to  another.  He  proposed  that  across  this  region  a  non- 
adiabatic  process  accurs,  rather  than  an  adiabatic  one, with  heat  being 
exchanged  among  the  particles  of  the  fluid,  but  none  being  received  from 
outside  it. 

Although  Rankine ! s  work  is  compatible  with  the  principle  of  conservation 
of  energy,  Rayleigh^  and  Hugoniot^  were  the  first  to  point  out  clearly 
that  inentropic  transition  in  a  shock  would  violate  this  principle. 

Hugoniet  showed  in  fact,  that  for  non-viscous  flow  without  heat  transfer 
(eutside  the  discontinuity)  entropy  must  be  conserved  for  continuous  flow, 
and  must  change  across  a  discontinuity.  Lord  Rayleigh  showed  that  dissi¬ 
pation  is  necessarily  present  in  a  shock  wave,  and  that  entropy  across 
auch  a  wave  must  hence  increase. 

5  B.  Riemann  -  Uber  die  Fortpflanzung  ^bener  Luftwellen  von  endlicher 
Schwingung3weite  -  Gesammelte  Werke,  1876,  p  144  or  Abhandlungen  der 
Ges ell shaft  der  Wissenschaften  zu  Gottingen.  Math-Phys  Klasse  8,  43, 
(1860) 

6  W.  J*  M.  Rankine  -  One  the  thermodynamic  theory  of  waves  of  finite  longi¬ 
tudinal  disturbance.  Phil -Trans.  A  160,  277-288  (1870) 

7  J.  Rayleigh  -  Aerial  plane  waves  of  finite  amplitude  Proc.  Hoy*  Soc.  A. 
84,  247-284  (1910) 

®  H.  Hugoniot  -  Sur  la  propagation  du  mouvement  dans  les  corps  et  3peciale- 
meat  dans  les  gaz  parfaits.  Journal  de  l’ecole  polytechnique,  5#,  1-125 
(1889)  Page  3 
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The  earliest  experimental  work  was  done  by  Mach^  and  hi3  coworkers,  who 
demonstrated  by  various  methods  that  shock  waves  spread  with  velocities 
greater  than  that  of  sound. 

The  following  discussion  is  a  detailed  theoretical  investigation  of  one- 
dimensional  shock  waves  and  derives  almost  all  of  the  results  obtained  by 
the  earlier  workers  in  the  field  by  application  of  the  principles  of  con¬ 
servation  of  mass,  momentum,  and  energy,  and  by  ingenious  mathematical 
procedures. 

The  investigation  of  one-dimensional  shock  waves  is  well  justified  in  a 
paper  which  is  supposed  to  be  primarily  concerned  with  spherical  shocks. 
Not  only  is  an  understanding  of  one-dimensional  shocks  necessary  for  com¬ 
prehension  of  spherical  shocks,  but  in  the  treatment  of  one-dimensional 
shocks  we  are  also  able  to  reach  some  conclusions  that  apply  to  shock 
waves  in  general. 

The  following  treatment  follows  that  of  Becker^  fairly  closely,  but  ad¬ 
ditional  mathematical  steps  and  interpretations  have  been  added,  and 
comments  by  Lewis  and  von  Elbe^“  have  been  included. 

9  e.  Mach,  Wiener  Berichte,  72(1875),  75(1877),  77(1878) 

10  R.  Becker  -  Stosswelle  und  Detonation.  Zeits.  f.  Phys.  8,  321(1921) 

11  " Combustion,  Flames,  and  explosion  of  Gases1'  B.  Lewis  and  G.  Von  Elbe, 
Academic  Press.  Inc.,  New  York,  1951 
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QUALITATIVE  DISCUSSION  OF  SHOCK  WAVE  FORMATION 
Before  proceeding  with  a  detailed  mathematical  analysis  of  shock  wave 
formation,  we  might  do  well  to  become  more  acquainted  with  that,  phenomenon 
on  a  moro  qualitativo  basis. 

Let  us  consider  a  long  tube  into  which  a  piston  is  inserted  at  the  left; 
the  tube  being  filled  uniformly  with  the  gas  of  the  surrounding  atmosphere. 
If  we  impart  &  small  velocity,  dw,  to  the  piston,  this  movement  causes  a 
weak  compression  wave  to  travel  to  the  right  with  the  velocity  of  sound. 
Then,  at  any  given  instant,  the  g&3  to  the  right  of  the  wavefront  is  un¬ 
changed  and  undisturbed,  while  the  gas  between  the  wave-front  and  the 
piston  is  compre33ed  an  amount  dp  and  has  the  velocity  dw.  If  we  now 
increase  the  velocity  of  the  piston  by  another  increment  dw,  a  second 
weak  compression  wave  will  proceed  to  the  right.  By  frequent  repetition 
of  this  procedure,  the  pi3ton  is  brought  to  its  final  velocity,  w. 

The  compression  waves  that  are  started  later  advance  with  greater  velocity 
since  the  velocity  of  sound  is  greater  at  higher  densities  (or  at  higher 
temperatures,  if  we  assume  the  compression  to  be  adiabatic),  and  since 
the  gas  itself  also  has  higher  flow  velocity.  Hence,  as  newer  waves  "catch 
up"  with  the  older  ones,  the  wave-front  becomes  increasingly  "steeper"  - 
when  older  and  newer  waves  are  completely  merged,  a  shock-wave  with  an 
extremely  large  pressure  gradient  is  formed. 
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A  MATHEMATICAL  INVESTIGATION  OF  SHOCK  WAVES 


Hie  Formation  of  Shock  Wavc3 

We  consider  a  tube  as  described  in  the  previous  section,  but  with  unit 
cross-sectional  area.  Let  denote  the  thickness  of  a  very  thin  disk 
of  the  gas  in  the  tube,  let  x  be  the  space  coordinate  measured  along 
the  tube,  and  let  t  denote  time.  Then  we  see  that  any  property  G 
of  a  material  particle  along  the  tube  is  a  function  of  x  and  t  only, 


Hence  we  may  write 

(1) 

dk-t  * 


'ICz 

*  3* 


where  u  =  denotes  velocity  along  the  tube  axis.  Further,  we  may 

show  (see  Appendix  X)  that 

cL-t  ^  3X 

Since  we  have  defined  K  as  very  small,  we  may  consider  all  properties 
constant  within  this  thickness.  If  p  denotes  density,  p  pressure, 

E  internal  energy  per  unit  mass,  we  note  that  £  contains  the  mass 
p  £  ,  the  momentum  up  ^  ,  and  the  energy  p  £  (  E  Letting 

Pt  represent  the  mean  pressure  acting  on  the  surfaces  of  our  disk,  and 
letting  X  be  a  heat  transfer  coefficient,  and  jjl  one  of  friction,  them 
from  the  principles  of  conservation  of  mass,  mementum  and  energy,  one 
may  write  the  following: 

(2->  Sr(ps'.-° 

(a» 

(2c)  u  *!£)] 


where  pt  =  p(  and  where  the*  relation  between  jll  and  the  ordinary 
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▼lscoaity  coefficient  'rp  is  17  -  \ /L  as  dictated  by  the  symmetric 

12 

property  of  the  pressure  tensor. 

The  above  three  relations  are  easily  modified  into 

(3»)  =  -  P  — 

w  J  dt  I  3x 


f5  (  P  A 


(  p  -ji  1  dg. 

'R  ^9x1  Px  dLt 


i  A  31 

p  3x  (  A  Dx 


The  manipulations  involved  are  shown  in  appendix  II.  If  we  introduce 


entropy  S  through  the  relation 

T  as  *  cdk  E  ^  p  edor  *  ol  E 
then  equation  (3c)  may  be  expressed  as 


P  dg. 

p  P- 


(4)  „  t  dS_  =  +  —/a  dl  \ 

p  +  c)x  l  A  sxJ 

which  expresses  the  variation  of  S  with  time,  under  the  influence  of 
yU.  and  X  ,  i.e.  friction  and  heat  conduction. 

Now,  let  us  neglect  heat  transfer  and  friction,  for  the  sake  of  discussion 
This  amounts  to  setting  A *  and  hence,  from  equation  (4)>  *=  0 

This,  in  turn  implies  that  2>  *  constant,  or  that  the  flow  occurs  adia- 
b&tically  and  reversibly. 


Applying  (1)  to  p  ,  and  using  (3a)  we  find 

<?a)  *  *  ?x  +  p  ~  ° 

Similarly,  applying  (l)  to  u.  ,  and  making  use  of  (3b)  after  setting 

/-o. 


12  Weber  and  Gaus,  Report,  d.  Phys.  I,  1,  p.  439 
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If  we  postulate  a  gas  for  which  p  =  p(f)  ,  (say  an  ideal  gas  for  which 
p  =  axpK  where  a  is  a  constant  and  k  the  ratio  of  heat  capacities,) 

the  above  becomes 
(5b) 


t)u 


ax 


_l_  2(1  ^  o 

p  ^P  ax 


Equations  (5a)  and  (5b)  may  be  solved  as  follows,  to  yield  solutions  of 
the  form  u(x^,p(*,t).  We  consider  the  X-t  plane,  in  which  the  direction 
of  any  line  element  (d'\*t)  may  be  given  by  *f  ,  where  f  is  defined 
by  ax*  if  at.  Along  this  line  element,  any  function  GM  changes  by 


the  amount 


ciC* 


( 


in 

ax  ' 


t-j 


’“•  v  H- ) " 


If  G  is  chosen  as  G-u^(f),  where  f  is  any  function  of  p  only,  with 
derivative  ^  3  we  find 

d  t  u  +  *  (  §x  f  +  i' 

By  adding  and  .attracting  ugj  ,  (V  within  th.  f*r.nth,„.a,  th. 

above  may  be  transformed  into 

dta+il  «  *  V 

Referring  to  (5a)  and  (5b),  we  find  that  the  right  side  of  this  expression 
vanishes  if 

u-w-h'-  p  ^  ~yui  hr  ■  ^ 

that  is  either  if  _  , - 

<’  ■  i  {Ifr  *~A  t  *  “■  *  (lap 

""  <'"p^  w 

We  now  define  for  ease  of  notation,  c*  =>  .  With  this  definition 

the  foregoing  results  may  be  written  more  concisely  as 

AVI  d  a  Ua  *  C 

But,  since  these  are  solutions  which  make  =  o  ,  we  find  that 

(64) 


tnese  are  aoiutions  muou  \  -  \  • 

p 

j  —  dp  =  Const,  along  curves  =  u. -v  c. 


and 

(6b) 


I 

j  dp  ^  Const  -  along  caf 


yes 


dt 


•c*  LL  -  CL 


?• 
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po  here  denotes  a  constant  of  integration,  which  is  evaluated  by  appli¬ 
cation  of  initial  condition^  and  found  to  be  the  density  of  the  undis- 
turbed  medium* 

If  we  introduce  Riemann's^  notation,  to  ijnprove  on  Becker^s  cumbersome 

system,  and  define 
p 

m  as  indicated  by  ft.  H.  Cole^*^ 

'  I#  ‘ 

the  above  equations  may  be  written  more  simply 
u-t  o-  =  Const,  along  curves  -  a.  ±  c. 

Now  we  shall  apply  these  results  to  the  piston  in  our  infinite  tube.  We 
assume  the  piston  at  rest  at  0  and  t -O  ,  let  it  be  accelerated  uni¬ 
formly  with  an  acceleration  a  for  a  time  T  and  finally  let  it  continue 
with  constant  velocity  a,  *  If  we  denote  quantities  pertaining  to  the 
piston  by  the  subscript  s  ,  the  motion  of  the  piston  will  be  defined  by 

C  «  k  ^  **  ,  ^4  =*  O  ts  T 

(7)  \ 

I  Xs  »■  a.Tfcs  +  i  *•  cx  T  *  -  T 


This  motion  is  shown  as  curve  Cl  in  the  X~t  plane  (Fig*  1)* 

In  the  entire  tube,  U=0  ,  p  *  po  at  time  ±*0  .  Since  u-fr  then  has  the 
same  value  everywhere  on  the  x-axis  (at  t»D  ),  and  since  the  curves 
(6b),  originating  from  the  positive  x-axis,  fill  the  entire  space  be¬ 
tween  C  and  the  x-axis,  then  in  that  entire  domain,  *  *  const* 
Similarly,  along  the  curves  (6a),  U.+  -  const.  Hence,  along  the 

curves  (6a)  both  u  and  CT  must  be  constant,  or  both  u.  and  p  must  be 

13  R.  H.  Cole  -  Underwater  Explosions*  Princeton  University  Press, 
Princeton,  New  Jersey,  1943 
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constant.  (  p  ~  const,  implies  =  const.)  But,  along  the  x-axis, 
(i.e.  att=D  ),  a=0  ,  p  *  p0  .  This  means  that  along  the  x-axis 
^  =  0  also,  and  hence  a-0~  =0  ,  so  that  the  relation  between 

u>  and  p  that  is  valid  everywhere  is 

(8)  M.  =  o-  =  [  P  j;  J|e  Af> 


Immediately  at  the  piston  (i.e.  along  the  curve  C),  is  given,  and 
hence  also  ps  by  the  above  equation.  Then  from  every  point  (  X*,  \  ) 


on  C,  the  line 


x-  x5  ^ 


(t  -ts)  (  U 


on  which  p  and  u.  have  the  constant  values  ps  and  ^may  be  drawn. 

For  piston  motion  given  by  (7),  the  portion  of  the  \-t  plane  between 
the  curve  C  and  the  x-axis  is  divided  into  three  regions  by  the  two 
lines  (9)  corresponding  to  1^-0  and  \  =  T  #  In  the  lowest  region 
u  *  O  ,  in  the  middle  (shaded)  region  u.  varies  from  O  to  u.l  ,  and 
in  the  upper  region  u.  =  ll,  (constant). 


In  a  gas,  where  p  =  <**  p  ,  =  o?  VC  pKM  9  whence 


v  p 

rp  pK-iM 

^  =  P  =  ^  \p. 

If  we  now  denote  by  c0  the  velocity  of  sound  in  the  gas  at  the  initial 

i — 

conditions,  c.0  =  <MK  p.  *■  ,  and  we  obtain  from  (8); 

!<x.  «  (t  -  c#)  ,  and  from  this 

t-  (£)*-(*)*-'♦  *£* 

The  slopes  of  curves  (6a)  and  (6b)  then  become 


U»  -v  C  *  u  -4-  u.  c, 


Co  ^T1  ^ 


M.  -c.  -  u.  — 


IA.  ~  <L0  = 
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It  is  important  to  notice  that  the  solutions  we  have  obtained  cease  to 
be  valid  as  soon  as  two  of  the  lines  (9)  intersect,  since  at  that  point 
of  intersection,  a  would  be  required  to  have  two  values  at  the  same 
time.  Hence  the  analytic  counterpart  of  one  wave-front  over-taking  the 
other  is  the  intersection  of  two  curves  (9).  The  place  X  ,  and  time 
T  of  this  over-taking  are  found  by  simultaneous  solution  of  (9)  with 

at  -  -  co 

which  is  obtained  by  differentiating  (9)  with  respect  to  ,  after 
making  U3e  of  (11)  and  (7)  for  replacing  ,  by  functions  of 


.  Then,  we  find 

,  X=  to T 

For  a  =  ZOO  "Ate.  ,  c„  =  330  "Aei  ,  K  =  i  4  ,  T=  O.s  vec  , 


find  '-‘•i  =  100  "}s<c  .  The  first  discontinuity  will  occur  at  X  «■  453m, 

T  =  1.3S  sec  (for  "t*-  O  ). 

Also,  by  (10),  p  p 

—  •*  i.s-t  SL  .»  1.34 

po  '  p- 

Fig.  2  corresponds  to  the  above  example.  It  represents  the  distribution 
of  velocity  along  the  axis  of  the  tube  at  the  times  indicated.  The 
"stiffening"  of  the  wave-front  is  easily  recognized  by  the  increasing 
steepness  of  these  curves.  Although  Becker  gave  a  curve  of  this  type. 
Fig*  2  is  plotted  from  my  computations.  I  plotted  various  lines  of 
the  type  (9),  making  use  of  (7)  and  (11),  in  Fig.  3.  From  this  figure, 
Fig.  2  is  obtained  easily  by  a  cross-plot. 


It  is  evident  from  the  foregoing  discussion,  that  the  occurance  of  a 
compression  shock  depends  on  the  condition  that  within  an  adiabatic 
wave-formation  the-  waves  traveling  in  a  denser  medium  tend  to  overtake 
those  in  a  less  dense  medium.  This  means  that  the  velocity  given  by 


(6a)  ^  ^  cxWF 

at  ^  Ha 


If  we  apply  (8),  the  above  condition  becomes 


must  increase  with  increasing  density, 
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Bat,  the  above  expression  may  be  written  also  as 

F  -  f 

or,  if  we  replace  p  by  ,  we  obtain  (as  demonstrated  in  appendix 


IV)  that  this  expression  may  be  written  as 

— r^~  .  -  —  ■  P  Si ^ ce  P-  <0  ,  <x*SL<i  \ y  O  4K^ 

j  <*ip  ava  *  a ir  J  J 

Condition  for  the  formation  of  compression  shocks  becomes 

(  >  O' 

V  <4  Tr  /alidb«*.4ic 

From  the  foregoing  we  draw  the  conclusion  that  in  a  given  medium  either 
only  compression  or  only  rarefaction  shocks  can  occur,  not  both,  depend¬ 
ing  on  whether  /i* LF  \  is  positive  or  negative. 


Necessity  of  Considering  Heat  Conduction  and  Friction 
As  mentioned  before,  the  equations  (5)  are  only  applicable  where  heat 
transfer  and  viscosity  may  be  neglected.  Hence,  the  foregoing  solutions 
obtained  from  these  equations  are  subject  to  the  same  limitations.  Since 
no  real  fluids,  are  entirely  devoid  of  these  effects,  equations  (5)  must 
yield  wrong  answers  when  temperature  gradients  (  )  or  rate  of  change 

of  volume  (  —  ~  vrx  )  exceed  certain  limits.  This  is  evident  from 

inspection  of  equations  (5).  These  magnitudes,  however,  tend  to  exceed 
any  finite  value  in  a  shock  wave,  according  to  the  foregoing  discussion. 
Hence,  use  of  equations  (5)  is  possible  only  until  the  discontinuity  occurs 
but  is  not  valid  at  the  shoe  k  wave . 

Intuitively,  when  the  wave  front  reaches  a  certain  steepness,  a  further 
increase  in  steepness  will  be  prevented  by  friction  and  heat  conduction; 
the  flattening  tendency  will  be  just  sufficient  to  counteract  the  tend¬ 
ency  for  steepening.  The  details  of  this  interaction  are  too  complicated 
for  mathematical  description,  since  various  sedondary  waves  are  formed  in 
the  process.  We  hence  must  be  satisfied  with  investigation  of  the  shock 
wave  after  it  has  become  quasi-steady. 
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After  the  shock  wave  has  reached  a  quasi-steady  state,  we  are  able  to 
investigate  it  using  a  coordinate  system  that  moves  with  the  shock.  Our 
problem  reduces  to  integration  of  equations  (3)  for  the  case  where  all 
derivatives  with  respect  to  time  vanish.  7/e  hence  replace  dt.  fay  *  -  x 
and  obtain  the  following  equations  from  (3)r 

“  +  p  =  c 


(>xX  ox  v  ^  p  ~  _  ° 

—  .  ,vf  - *  *  3*  c  A  &) 

One  may  integrate  these  equations,  as  I  have  shown  in  appendix  V,  to 
obtain  the  following  results. 

(12a)  u  «  r MV* 

.  .  ...  ^  ar. 


M  ir  -r  p 


-  J  K  M 


p  J  ir  —  -V  M  ^  <r  ^  -  P 


a  a  t_ 

FA  ■  cA  } c 


where  M,  J,  F,  are  constants  of  integration,  and  energy  E  and  tempera¬ 
ture  T  may  be  considered  as  functions  of  f  and  ^  . 

We  seek  relations  between  the  values  ,  f.*,  .  q  ,  before  the  shock  front 
and  the  values  p.,  ,  behind  it.  This  is  easily  accomplished -when 

we  realize  that  only  within  the  wave  front  do  and  differ  markedly 
from  zero.  Hence,  for  any  place  outside  the  wave-front 

Lp  . 

(13)  MV+  p  -  iF  *  p  '  J 

E  -v  ^  **  P  ~ 

Comparing  two  such  places,  denoted  by  subscripts  1  and  2, 

_  Uc 

(14a)  irq  ~  ir^_ 

(Hb)  +  p.  -  V  +  Pi 

wl 

(14c)  C,  ■*  ^4*  +  F,  iTi  -  1=2.  +  ^  +  Pi  irz. 


Up 

*  Pi 

Pi 

e, 

-t  — ^ 

+  p,  «T,  - 

£2. 
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These  fundamental  equations  are  hence  independent  of  friction  /i  and 
heat  transfer  They  are  identical  to  relations  obtained  by  requir¬ 

ing  conservation  of  mass -t  momentum,  and  energy,  across  the  shock  wave.^ 
Hence,  the  mathematical  artifice  of  introducing  a  discontinuity  is 
justified*  (However,  by  this  means  we  obtain  no  insight  into  the 
processes  within  the  wave  front). 


Relations  obtained  for  Compression  Shocks 

In  order  to  obtain  a  clearer  view  of  the  meaning  of  equations  (14),  we 
solve  (14a)  and  (14b)  for  ,  and  ,  then  substitute  the  results  into 
(14c),  We  obtain 


(15a) 


*  -  F\ 


(15b) 


^  ~  TT. 


fV  -  P, 


(15c)  Ez-t,  =>  -fe  (p,  +  pj(  v,-n) 

Equation  (l5c)  is  known  as  the  Hugoniot  equation.  Note  that  this  is 

^  , 
different  from  =  I  pdv~  y  as  would  pertain  to  isentropic  flow,  but 

>vx 

that  for  small  enough  differences  of  £  and  vr  the  above  approaches 

Z  p4\Tf  which  is  the  isentropic  relation*  It  may  be  shown^  that 
for  the  same  volume  change  (15c)  will  read  to  a  larger  =  than 
the  corresponding  isentropic  relation. 

The  velocity  of  propagation,  U  ,  of  the  shock  wave  into  a  medium  at  rest, 
and  W  ,  the  velocity  of  flow  behind  the  wave  front  are 
(16)  C  LJ  =  u,  =  xr, 

<  I  ..  _ _ _ 

(  W  *  U,  -U.  <.  =  t  V,  -  iro  J  ~ 

Hence,  given  the  state  of  the  undisturbed  medium  (  ),  and  the 

"shock  pressure",  p^,  plus  an  appropriate  equation  of  state,  we  are 
able  to  compute  all  other  values. 


El  -  E 


14  This  is  done  by  Lewis  and  Von  Elbe,  op.  cit* 
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We  illustrate  this  for  an  ideal  gas,  where 


(17a) 

p  TT  ■  * 

■RT 

and 

(17b) 

E, 

-  (  U  -T,  ) 

between 

X.  and 

T,  . 

Vie  let  _ 

(18)  P'  -  >7- 

(  TT  -  p*/p. 

and  obtain  (Appendix  VI) 

(19a) 

TV 

TT+ 

TT  — - - 

TT  ' 

(19b) 

y% 

TT  -*•  f 

f.  TT  -v 

(19c) 

T>^  - 

TT  -  £  ,  *  » 

S,  -1 

(19d) 

W* 

*  p.-w,  (  fc;  -  »  ) 

where  cv  is  the  mean  heat  capacity 


(19d)  W*  -  P,^,  <  t,  ~ »  )  -"T*  - 

For  a  perfect  gasl5,  c*  =  hence  ^  •  For  a  diatomic  perfect 

gas15,  <  =  \-4  ,  and  hence  6  .  If  TT  is  large  compared  to  6,  (19a) 

indicates  that  T  increases  proportionally  to  p  .  Hence,  we  must  intro¬ 
duce  $  ,  as  a  function  of  T  . 


Becker  calculated  the  figures  given  in  Table  1,  using  .  r  4  -v 

which  he  supposed  to  be  good  up  to  about  3000° C.  (This  exp~< 

sion'^valid  f or  *Ot  and  COt ,  was  taken  from  Pier,  Zeits.  f.  Elektrochem. 
15,  536  (1909)  and  16,  897  (1910),  and  Siegel,  Zeits.  f.  phys .  Chem.  87, 
641  (I914)i)  Although  Lewis  and  Von  Elbe  state  that  this  expression 
is  somewhat  in  error,  Becker* s  work  still  serves  to  point  out  relative 
magnitudes. 


15  Obert.  Thermodynamics,  Mc-Graw-Hill  Book  Company,  Inc 

New  York,  1946.  p.  167,  equ.  (7-12) 

Page  15 


301  017 


and  with 


With  the  above  we  obtain  «  5 .82  +  O.46X40'1 
this  ^  ,  (19a)  becomes  a  quadratic  in  Tt  #  The  other  relations  (19) 
then  yield  all  other  required  values,  X  =  0°C. 


TABLE  I 


TT 

TV 

T. 

V, 

TV 

U 

0%^ 

W 

i 

(Vi 

2 

1.23 

1.63 

336 

452 

175 

1.63 

330 

5 

1.76 

2.84 

482 

698 

452 

11.4 

426  • 

10 

2.56 

3.88 

705 

978 

725 

34.9 

515 

50 

6.28 

6.04 

2260 

t 

2150 

1795 

296 

794 

100 

14.15 

7.66 

3,860 

3,020 

2590 

I 

699 

950 

500 

44.8 

11.15 

12200 

1 

6570 

5,980 

5560 

1 

1433 

1000 

70.0 

14.3 

19100 

9,210 

8560 

14300 

1 

3,710 

2000 

106.2 

18.8 

25j000 

12900 

12210 

37,600 

2p70 

3000 

134.4 

22.3 

36700 

15,750 

15050 

1 

66900 

2180 

The  total  impulse  \  , 

which 

appears  in  the  next-to-last 

column,  is 

an  important  measure  for  evaluating  the  effect  of  a  shock  wave  on  an 

obstacle*  t  is  composed  of  the  static  pressure  difference  p*.  -  P» 

"2. 

and  of  the  momentum  of  the  gas  behind  the  wave-front,  p*.  w  •  With 
the  value  of  VV  from  (16),  and  by  substituting  =  ^  , 

i  ~  -  (>zW2»  tFx-f.  )  w  A  .  (TT-.) 

In  order  to  illustrate  the  foregoing  calculations  for  a  liquid,  we  use 
Tammann's  equation  of  state^,  which  is  valid  up  to  very  high  pressures. 

16  G.  Tanmann  -  Ann.  Physik  37,  975  (1912) 

Page  16 

3°i  Q.i& 


It  is 


(20)  p  -  ^  * 

where  c,  b,  and  K,  are  constants.  The  general  relation  for  internal 
energy,  dE.  =  Cv  dT -  P)  av  then  may  be  integrated 
to  yield 

(20a)  E  =  T  + 

If  we  define 

(21)  p'  =.  p  +  K  '  v-'  =  IT  -  to 

and  apply  (20  )  to  the  Hugoniot  equation  (15c),  we  obtain 
C-ld  -  T,  )  -  t  (p,'’pi)(v/  -  V'z) 
and  (20)  is  transformed  into 

p'  \r*  -  C.T*  a 

The  above  equations  are  structurally  identical  to  those  for  a  gas,  so 


that  the  solutions  (19)  can  be  carried  over  directly*  Here, 

*v=  7~:  -V  \  ,  and  IT'  ='  ^  ,  making 

’  °  P,  +  K 


(22)  S.  ,  TV 


rr\+\ 


-to 

•  Vt  -to 


-n-1^  <, 


V,  -b 


Table  II  applies  to  ethyl  ether,  for  which  K-  2129  .  c-~  •,00' 

b  -  94  ~  \  v,  =  t.  It  ^  =•  .  S64  ,  p,  =  l  ,  T  -  c 


T/sbi-s  TE. 


p*.  C  oJrwv) 

®C"V4«J 

10G 

1.6 

1.6 

1260 

1,000 

15.6 

15.6 

1445 

10,000 

85 

113 

2680 

20,000 

123 

211 

3000 

60,000 

201 

594 

5010 

100,000 

245 

975 

6430 
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Structure  of  the  Shock  Wave 

Becker  investigated  the  structure  of  the  wave-front  by  carrying  out 
the  integration  of  equations  (12)  for  gases  in  which  **  and  /c  ful¬ 
fill  certain  requirements*  The  usefulness  of  his  calculations,  of 
course,  depends  on  how  well  a  gas  under  consideration  meets  these 
requirements.  For  air,  Becker  admitted  an  error  of  about  1($  in  his 
calculations  for  v  ,  ,  Ti  in  the  shock  front,  but  he  computed  thick¬ 

nesses  of  shock  fronts  on  the  order  of  magnitude  of  molecular  mean  free 
paths*  However,  for  such  magnitudes,  the  physics  of  continua  (hydrodynami 
can  not  apply,  and  Becker  admitted  his  theory  was  inadequate* 

Becker  obtained  the  following  shock  front  thicknesses  L  ,  for  air  with 

=  1  atm,  m  0°C 

tatra  )  2  5  10  100  1000  2000  3000 

447  117  66  16.5  5.2  3.6  2.9 

(According  to  Kinetic  theory  of  gases,  the  mean  free  path  at  1  atm, 

0°C  is  about  90  X  1CT^  cm,  and  the  mean  distance  between  molecules  is 
3.3  X  10"7  cm). 

But  Becker's  theory  was  vindicated  by  L.  H.  Thomas-*-7  in  1944.  He  showed 
that  all  shock  waves  in  air  are  a  few  mean  free  paths  thick,  by  consider¬ 
ing  the  increase  in  the  coefficients  of  thermal  conductivity  and  vis¬ 
cosity  with  increasing  temperature  and  pressure,  which  Becker  neglected, 
and  by  application  of  Kinetic  theory. 

Experimental  substantiation  of  theories  of  processes  within  a  shock  front 
Is  practically  nil.  The  extremely  transient  nature  of  shock  front  passage 

^  L.  H.  Thomas  -  Note  on  Becker’s  Theory  of  the  Shock  Front.  J*  Chem. 
Phys.,  12,  pp  449-57  (November  1944) 
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makes  instrumentation  a  formidable  problem,  and  besides  makes  knowledge 
of  conditions  within  the  shock  front  of  little  practical  significance. 

G.  R.  Cowan  and  D.  F.  Hornig13  developed  equations  for  computation  of  the 

density  profile  of  a  shock  wave  from  measurements  of  reflectivity.  They 

py  _ 

also  performed  some  measurements  in  nitrogen,  keeping  a  constant  P, 
1.71,  and  varying  p*  .  The  following  are  their  results: 

Pi  5.78  4.63  2.86 

i  oo-\^y  iso  200  320 

The  accuracy  of  their  measurements  is  estimated  at  25%;  hence  the  above 
values  are  slightly  larger  than  those  calculated  by  the  Becker  -  Thomas 
Theory.  Thi3  discrepancy  may  be  due  to  the  assumption  of  hard-sphere 
molecules  made  by  Thomas,-  which  is  not  too  good  for  a  diatomic  gas  like 
nitrogen.  But  assumptions  made  by  Cowan  and  Hornig  may  possibly  lead  to 
&  greater  error  in  their  work  than  was  estimated. 


G.  R.  Cowan  and  D.  F.  Hornig,  The  Experimental  Determination 
of  the  Thickness  of  a  Shock  front  in  a  Gas. 

J.  Chem.  Phys.  18,  1008-18  (August  1950) 
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Section  II 

THE  PROPAGATION  OF  SPHERICAL  SHOCK  WAVES 
INTRODUCTION 

Although  Lord  Rayleigh^  solved  the  hydrodynamic  equations  for  plane 
shocks  analytically  (assuming  adiabatic  relations  between  pressure 
and  density),  such  an  approach  can  not  be  taken  in  the  case  of  spherical 
shock  waves.  The  difficulty  arises  from  the  spherical  divergence  terms 
that  appear  in  the  partial  differential  equations  that  describe  spherical 
notion. 

A  solution  of  the  problem  may  of  course  be  obtained  by  numerical  integration 
of  the  pertinent  partial  differential  equations.  Penney2^,  and  later  Penney 
and  Dasgupta2^  hAve  carried  out  calculations  of  this  kind  for  spherical  TNT 
charges.  They  employed  a  method  of  integration  based  on  the  equations  of 
Riamann^,  which  greatly  simplifies  the  numerical  work  and  reduces  the  re¬ 
quired  labor..  But,  even  with  this  simplification,  this  method  is  extremely 
tedious  and  complex.  Furthermore,  the  Riemann  equations  apply  only  where 
dissipation  effects  can  be  neglected.  The  desirability  of  developing  more 
rapid  and  flexible  methods  is  hence  evident, 

J.  G.  Kirkwood  and  H.  A.  Bathe22  developed  one  such  method,  applicable  to 
underwater  explosions  of  spherical  charges.  However,  they  used  the  approxi¬ 
mation  ef  adiabatic  flow,  neglecting  the  increase  in  entropy  at  the  shock 
front.  This  assumption  is  warranted  for  water,  because  the  entropy  in¬ 
crement  produced  by  &  shock  wave  in  liquids  is  very  small.  However,  for 
gases  the  entropy  increase  is  not  neglible. 

^  W#  G.  Penney,  British  Report  RC  -  142  (1941) 

^  W.  G.  Penney  and  H.  K.  D&sgupta,  British  Report  RC  -  333  (1942) 

22  J.  G.  Kirkwood  and  H,  A.  Bethe.  -  Pressure  waves  produced  by  an  under¬ 
water  explosion,  I,  OSRD  5SS  (1942) 
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Another  theory  of  underwater  shock  waves  has  been  presented  by  Osborne 
23 

and  Taylor  J  •  This  theory,  however,  is  based  upon  the  acoustic  approxi¬ 
mation,  and  is  hence  strictly  valid  only  for  small  excess  pressures  at 
relatively  large  distances  from  the  source.  This  approximation  is  also 
not  valid  for  very  large  distances  from  the  source,  as  indicated  by 
G.  B.  Whitham2^. 

The  following  theory,  which  we  shall  examine  in  detail,  was  developed 
by  S.  R.  Brinkley,  Jr.,  and  J.  G.  Kirkwood2^.  it  is  superior  to  those 
discussed  above,  and  to  others  developed  for  shock  waves  in  air,  which 
have  been  based  on  ideal  gas  adiabatics  with  constant  heat  capacity^.  . 
Brinkley  and  Kirkwood  account  in  their  theory  for  the  finite  entropy 
increment  occuring  in  the  fluid  due  to  passage  of  the  shock  wave.  Hence, 
this  theory  is  equally  valid  for  the  propagation  of  shock  waves  in  liquids 
and  in  gases. 

The  development  which  follows  is  essentially  that  of  Brinkley  and  Kirk¬ 
wood;  however,  I  have  specialized  it  for  application  to  spherical  shocks 
(since  their  discussion  applies  to  plane  and  cylindrical.flow  also) and 
I  have  inserted  additional  explanation  and  intermediate  mathematical 
steps  where  these  seemed  applicable. 

23 

F.  M.  Osborne  and  A.  H.  Taylor.  -  Nonlinear  propagation  of  underwater 
shock  waves.  Phys.  Rev.  70,  322  -  8  (September  1,  15,  1946) 

G.  B.  Whitham.  -  The  propagation  of  spharical  blast.  Proc.  Roy.  Soc. 

A  ,203,  571  -  81  (October  24,  1950) 

25 

S.  R.  Brinkley,  Jr.,  and  J.  G.  Kirkwood.  -  Theory  of  the  Propagation 
of  Shock  Waves.  Phys.  Rev.  71,  606  -  11  (May  1,  1947) 
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THE  BRINKLEY  AND  KIRKWOOD  THEORY 


Brinkley  and  Kirkwood  use  a  hybrid  form  of  the  hydrodynamic  equations, 
utilizing  the  Euler  form  of  these  equations,  but  introducing  Lagrangian 
Coordinates^.  In  order  to  derive  the  equations  as  they  are  used  in  this 
theory,  we  define  the  following  symbols. 

Let  r  be  the  Euler  coordinate  at  time  "t  of  an  element  of  fluid  with 
Lagrange  coordinate  .  Let  p  and  p  denote  pressure  and  density, 
and  let  the  subscript  (0)  refer  to  the  undisturbed  fluid.  The  Euler 
velocity  of  sound  c  ,  is  given  by  cJ*  s  (  ,  and  let  P*  =  p  -  p>  . 


The  equation  of  motion  (conservation  of  momentum)  in  the  Euler  form,  may 
be  written  in  vector  notatiorTas 

(23)  P  <^-  +  7'p  *  0  ,  where  \*  is  the  velocity  vector 

For  spherical  flow,  all  quantities  are  functions  anly  of  radius  r  and 
time  t:  ,  and  the  velocity  is  purely  radial.  For  this  type  of  flow,  (J23) 


becomes 


<»> 

where  is  the  particle  velocity,  or  ■  u.  = 

But  cl  =  ( R,*fc)  and  Euler’s  notation 


,  and  Euler’s  notation  means  R  *  const. 


Hence 

diu  (  \  /  <)^.  \  I  Ou.  \  f 

dL*b  =  )  l  'D-b  V  t>-t  l-pi.  *  \  )*£. 

We  note  that  a  spherical  shell  with  radius  R  and  thickness  dR  contains  a 
mass  of  fluid  given  by  4TTp0  d'R-  ,  (R  is  by  definition  the 
position  in  the  undisturbed  fluid  of  a  volume  element  which  at  time  t  is 


at  r  ),  and  that  the  same  mss  of  fluid  is  contained  in  a  spherical  shell 


26 

G.  I.  Taylor  The  air  wave  surrounding  an  expanding  sphere. Proc .  Roy. 

Soc  A.  186,  273  -  92  (September  24,  1946)  -  This  assumption  fails  badly 
near  the  explosive  source. 

27 

See  the  "Glossary  of  Hydrodynamic  Equations”  near  the  end  of  this  paper 
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of  radius  r  and  thickness  dr  y  which  contains  a  mass  4r 


Hance,  at  time  h  $ 


which  wo  find 


Then,  considering  p 


:  4-nrir'dr 

P.'R* 

p  r1 

p  («,*) 


,  from 


V  T>*7 


l£\  />t 


a*.)  £_e 

c)R/t  Po'R 


Substituting  the  foregoing  values  into  (24),  we  obtain 

■  Flg  .*  Ji  jlE.)  .o' 

or,  as  in  the  paper  by  Brinkley  and  Kirkwood, 

(25a)  -5!:  (  Oi)  ^  JL-fV?) 

r1-  Ut  /r.  p0  ~  0 

since  d"P  =  dp 

Euler's  equation  of  continuity  in  vector  notation  appears  as 
(26)  +-  p  (  V-  W-) '  =  0 

Again,  constant  R  is-  implied.  Considering  p  s  p  (r,  t)  again, 
then 

a-t  dp  d.t  O  d-b  c*L V  DK;t\  (.2>-fc)Ri  <-x\Z>*  J\ 

Further,  for  spherical  flow,  and  considering  u  =  u  (r.t), 

V/  =  ^  ^rCr^u.) 

by  the  definition  of  v  ,  for  curvilinear  coordinates2®;  or 
V-  v/  =  ( ]t±.\  .  hL.  _  /  3m.  \  /  OR  \  a.'M. 

^4  r  -  15^4 ( yr)* +  - 

=  _e£  / 

_ J’o'R-*  \ 

28 

See  for  instance,  H.  Lass,  "Vector  and  Tensor  Analysis"  McGraw-Hill 
Book  Co.,  Inc.,  New  York,  1950.  p  54,  equ.  (89) 
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Substitution  of  the  foregoing  results  into  (26),  and  division  by  p 

yields  Kirkwood  and  Brinkley's  form  of  the  equation  of  continuity: 
(25b)  .  , 

Summarizing,  we  have  the  new  equation  of  motion 

(25*)  ( 1>±\  j_  (  yp_\  _  ^ 

r'K-Z-tU  Po  \  2^4  ~  0 

and  the  equation  of  continuity, 

<25b)  p  (  lf)R  *  ^  .  o 

Eq^tlon.  (25)  «r.  to  6.  soiv.d  subject  to  th.  lnltlia  Saialtlon,  5p.c 

ri.d  «  .  curve  to  th.  R,  t-pl.„.  and  u  th.  Hugoniot  conditioh,  .t  th. 

shock  front  (Appendix  VII): 

—  p>«  u.  U 

(26)  P  (  U  -u)  -  p0U 

“  "  h*r'  dTOt*»  the  *"thaipy  increment  .gp.ri.nc.d  by  th. 

th.  fluid  in  trav.raing  th.  shock  front,  »nd  U  is  th.  velocity  of  th. 
shock  front.  (Th.  Hugoniot  conditions  ar.  ....  to  .pply  to  spherical 

•hock  fronts  Use,  sine,  i„  th.ir  d.riv.tion  no  assumption  ..  to  th. 
shap«  of  the  fronts  was  made). 

At  th.  shock-front  itself,  r  =  R,  so  that  equation,  (25)  for  the  shock 
front  simplify  to 


(~)  +  i.  1 

1  OPN 

'  <5  t  'ft.  p0  ' 

Or  J 

(25d)  ^  _i_  /OpN  2^ 

P°  V  3R  )t  ^  V  2 *  )r.  ^  ° 

We  define  now  a  derivative  in  which  the  shock  front  is  stationary  by  tl 
•perator  ^  ^  ^  ~  - 

^  ^  +  -Sr.:  5T  -  bh 

*nd  apply  this  operator  to  the  first  of  equations  (26). 

Page  24 


301  02$ 


After  transformation  we  obtain  (Appendix  VIII) 

(27) 

W  ^  po  po  U  2>  t  .  . 


where 


Z.  au 

O  ciP 


We  note  that  all  of  the  coefficients  of  equations  (25c),  (25d),  and  (27) 
may  be  expressed  as  functions  of  pressure  alone  by  means  of  the  Hugoniot 
conditions  (26)  and  an  appropriate  equation  of  state  of  the  fluid.  If  we 
dan  find  a  fourth  relation  between  the  four  partial  derivatives,  we  would 
be  able  to  solve  for  each  of  the  four  derivatives  as  a  function  of  P  and 
R  and  to  formulate  by  means  of  the  operator  we  defined  previously,  an 
ordinary  differential  equation 

(J8)  ^  •  ||  *  t|j  If-  >  *0 

dK  U  ox 

for  the  peak  pressure  P  of  the  shock  wave  as  a  function  of  the  distance 
R.  (Henceforth  we  shall  use  and  P*  to  denote  particle  velocity  and 
excess  pressure  behind  the  shock  front,  and  reserve  the  unprimed  letters 
for  these  quantities  at  the  shock  front). 

The  physical  basis  for  establishing  the  desired  fourth  relation  is  the 
fact  that  as  a  shock  wave  passes  through  a  fluid,  it  leaves  in  its  path 
a  residual  internal  energy  increment  in  each  element  of  the  fluid.  This 
energy  increment  is  determined  by  an  entropy  change  produced  by  passage  of 
the  shock  front.  As  a  result,  the  energy  propagated  by  the  shock  wave 
decreases  with  distance  from  the  source. 

Clearly,  the  work  done  by  the  shock  wave  after  it  passes  a  point  R  is 
ultimately  transformed  into  internal  energy.  Consider  the  work  done 
by  the  source  of  initial  radius  a.  „  This  work  may  be  resolved  into  two 
parts:  that  which  is  transformed  into  internal  energy  of  the  fluid  at 
a  pressure  ^  within  a  sphere  of  radius  R,  and  that  work  which  is 


Note  that  passage  of  a  shock  front  through  a  fluid  element  leaves  it 
at  increased  entropy  and  internal  energy,  from  which  it  returns  to  the 
initial  pressure  .  Absence  of  heat  transfer  is  assumed  throughout. 
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done  on  this  spherical  surface  by  the  external  fluid. 

The  first  part  of  this  work,  Wot  >  may  be  found  by  integrating  Ed m  be¬ 
tween  the  initial  and  final  spherical  surfaces.  The  mass  dm  contained  in 
a  spherical  shell  of  radius  r0  ,  and  thickness  ,  is  dm-  =r  4Trp<, 

i3  used  here  because  the  fluid  is  initially  in  the  undisturbed  state. 


Hence, 


W„,  -  $ 


t  dm  »  ^TT 


j  er< 


r?  <dr. 


This  integral  may 


be  evaluated  if  E  =  E  (  ),  is  known;  ~P  being  the  peak  pressure,  and 
if  P  =  (  f0  )  can  be  found. 

The  second  part  of  the  work,  Wb2_  ,  may  be  found  from  the  definition  of 
mechanical  work,  as  the  integral  of  normal  force  times  distance.  Here,, 
our  force  is  pressure  times  area  of  a  spherical  surface,  and 


f  °°  30 

W0i  -  4ir  r‘(p'+.ft) 

Hence, 

(29)  -R.  c- 

~  v  (  r'  u.1  at 

Ja-.  H.OO 

The  term  involving  p0  in  the  time  integral  gives  the  product  of  p\>  and 
the  volume  displacement  of  the  fluid  element  initially  at  R.  This  dis¬ 
placement  is  obviously  the  sum  of  the  outward  volume  displacement  of 

the  inner  boundary  of  the  fluid  (i.e.,  the  expansion  of  the  explosion  pro¬ 
ducts)  and  the  displacement  of  the  volume  of  fluid  initially  between  shells 
of  radii  ac  and  R  to  shells  of  radii  cl'  and  R*  .  We  may  hence  write 

-4TT  f  ft  r^u'ieUt  -  f0  aV3  +•  [.  (  ~  (  rgdv;  1 


t„(R) 


Note  that  p  -  P  /  jo  ,  and  udt  is  the  displacement  of  a  fluid  element 
in  time  dt  \  t#(R)  is  the  time  of  arrival  of  the  shock  front  at  R. 
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or 


*  o# 

4TT  I  p.Wdt  =  p.  *V3  *  4TTp.  (  (-&■-»)  d^. 

H.tR)  •’a.  ( 

from  the  relation  pr2dr  =  p0  r„2  dre  between  Euler  and  Lagrange 

coordinates  as  previously  derived. 


If  we  now  introduce  h  “  E  /  pQ  £*(-^),  the  specific  enthalpy  incre¬ 


ment,  we  find 


\  E  Po  r*  <drQ  a  [  Kp0  ro'diTo  -  (  f°  Po  ("3  ^ 


■ 


^  Po  — 


Substituting  then,  (29)  becomes 

Po  ^ 


j  ^  C^-o  ^  d'c« 

^<3.0 


(30)  Wo_ 

^nr 


4o 


a-Tp 


r'TV 

JU-*> 


dit 


The  time  integral  is  assumed  to  vanish  as  the  shock  front  passes  to 
infinite  distance.  Subtracting  from  (30)  its  value  at  R  =  ,  we 

obtain 


O  = 

from  which 


»  1  -  (  ^p0r0^dr0  +  (  ^ 


at 


(31)  (  Wo0fQxJ<0  -=.  (  r^'P’cJ  tjL-fc-  =  C^-) 

.  R  VO 

The  energy  of  the  shock  wave  at  the  point  R  is  by  definition  the  work 
done  on  the  fluid  exterior  to  R. (Since  total  energy  of 'shock  wave  is 
energy  expended  from  source  to  R  plus  energy  expended  outside  R,  the 
energy  remaining  in  the  shock  wave  when  it  is  at  R  is  the  total  energy 
minus  that  expended  between  the  source  and  R,  or  the  energy  expended 
outside  R).  The  shock  wave  energy  at  R  is  hence  given  by  4TT  JDC^)» 


D(R)  as  written  above  does,  however,  depend  on  the  size  of  the  source,  and 
it  is  hence  desirable  to  normalize  the  tine  integral  to  a  value  which  is 
independent  of  this  factor.  This  is  done  by  expressing  the  integrand  as 
a  fraction  of  ^its  initial  value  R2Pu,  and  choosing  a  reduced  time  scale 
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for  which  the  initial  slope  of  the  integrand  has  unit  value.  Inasmuch 

2 

as  the  observed  energy-time  curves  are  approximately  of  the  form  R  Pue  , 
it  is  convenient  to  use  as  a  time  unit  the  initial  logarithmic  slope  of 
the  curve,  defined  by 


k  (-H  pVr7t  =t; 


After  the  differentiation  is  carried  out,  this  becomes 


(32)  ^  _  Jl 

/jl  P  Dt:  ^ 

Then,  D(R)  may  be  expressed  as 


u_  t 


(33)  XXLR.3  -  -R’-  P>  -L> 


Where  ~V  -  ~u 


(r)  =  j  f  Cr, 


t)  ar 


T  being  given  by 


-t  -t0c^) 


\  CR.-rV  =. 


r  ^  P’u.' 
■R.  P 


Elimination  of  jx.  between  (32)  and  (33)  yields  the  desired  fourth  re¬ 
lation  between  the  partial  derivatives  at  the  shock  front: 


P  2>t 


i  Du. 
u.  3t 


The  set  of  equations  (25c),  (25d),  (27),  and  (34)is  exact,  involving 
integrals  of  equations  (25a)  and  (25b)  for  knowledge  of  the  reduced 
energy-time  function  f(R,'T). 

However,  if  f(o^y)  is  initially  a  monotone  decreasing  function  of  T, 
f(R,T)  will  remain  so,  aad  in  fact  will  at  larre  R  become  asymptoti¬ 
cally  a  quadratic  function  of  T  corresponding  to  the  linear  form  of  the 
pressure-time  curve  that  has  been  shown  to  be  stable  at  large  distances31* 


J.  G.  Kirkwood  and  H.  A.  Bethe.  03RD  588(1942) 
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This  means  that  -V  is  a  very  slowly  varying  function  of  R,  for  which 
sufficiently  accurate  estimates  for  many  purposes  may  be  made  without 
explicit  integration  of  the  hydrodynamic  equations.  The  assignment  of 
a  value  which  is  independent  of  R  to  v  is  equivalent  to  imposing  a 
similarity  restraint  on  the  energy- time  curve.  This  type  of  restraint 
has  been  used  before,  and  has  proven  to  be  a  very  reasonable  assump¬ 
tion.26^ 

Ibe  initial  energy-time  and  pressure-time  curves  of  an  explosion  wave 
are  rapidly  decreasing.  For  an  initial  estimate  of  ^  ,  an  expansion 
of  the  logarithm  of  the  function  in  a  Taylor  series  in  time  is  appro¬ 
priate,  (The  well-known  -peak-approximation).  This  corresponds  to  an 
exponential  ^l*T)  $  which  gives  the  result  v  -  j  T  =  I 

o 

For  the  asymptotic  quadratic  energy-time  curve,  tT)  =  (  >  - 

W  T  ~  2  o^d.  V  t~r )  =.  o  ipr<  ~r  >  z  , 

v  -  *■¥. 
o 

As  a  convenient  empirical 

interpolation  formula  between  these  two  values,  Kirkwood  and  Brinkley 
have  used,  in  a  different  paper, ^ 

-y  |  -  e  ~ 

in  a  series  of  calculations  of  the  peak  pressure  distance  curves  for 
spherical  shock  waves. 

Equations  (25c),  (25d),  (27),  and  (34)  may  be  solved  simultaneously  for 
the  four  partial  derivatives,  and  an  ordinary  differential  equation  for 
^  33  a  function  of  R  may  be  formulated  with  the  aid  of  equ.  (28)  . 


G.  Taylor  -  Similarity  solutions  to  problems  involving  gas  flow  and 
shockwaves.  Proc.  Roy.  Soc.  A.  204,  p  8-9,  (November  ly>0) 

33 

J.  G.  Kirkwood  and  S.  R.  Brinkley,  Jr.,  03RD  4814  (1945) 
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A  second  ordinary  differential  equation  relating  D  to  R  may  be  obtained 
by  differentiation  of  (31). 

The  results  of  the  above  manipulations  (as  indicated  in  appendix  X) ,  may 
be  written  in  the  following  form. 


(35) 

dP 

d-R 

-  ^  L(R) 

dP 

-v  ^  M 

dR  = 

where 

L  (V  ) 

=  po  W.  C 

M  C?) 

'  ^ 

PoU1-  2.0  + 9)- 

N  CP) 

=  A  (  Po/p)  f  2  C 

z  (  9)  _  ^ 

p  dO 


3  ~  1  cj  <4P 

The  functions  L(  P  ),  M(  V  ),  N(  P  )  can  be  evaluated  as  functions 
of  the.  pressure  by  means  of  the  equation  of  state  of  the  fluid  and  the 
Hugoniot  relations  (26),  (This  is  illustrated  for  the  case  of  a  perfect 
gas, in  appendix  XI),  The  equations  (35)  may  then  be  integrated  numerically, 
by  the  use  of  standard  methods^,  provided  the  initial  conditions  are  known. 

The  constants  of  integration  may  be  determined  in  one  of  two  ways: 

(1)  by  the  initial  conditions  at  the  boundary  of  the  explosion  pro- 
ducts  following  detonation,  or 

(2)  by  an  experimental  pressure-time  curve  at  a  selected  distance. 

The  first  method  is  a  theoretical  one,  and  is  based  on  the  thermo-chemistry 
connecting  the  explosive  and  its  products.  The  study  of  detonation  is  a  field 


See,  for  instance,  J,  B.  Scarbrough,  "Numerical  Mathematical  Analysis.rr 
Johns  Hopkins  Press,  Baltimore,  Md.,  1930,  pp218  ff. 
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in  itself,  and  is  hence  beyond  the  slope  of  this  paper.  3ome  discussion 
of  this  may  be  found  in  R.  H.  Cole’s  book  ’Underwater  Explosions”,  and 
Brinkley  and  Kirkwood  have  published  a  paper^giving  tables  and  graphs 
for  application  of  this  type  of  boundary  condition  to  numerical  integration. 
According  to  Cole,  the  most  convenient  initial  constants  in  this  case  are 
the  initial  pressure  and  either  the  parameter  Jl  or  the  total  shock-wave 
energy. 

If  the  semi-empirical  approach  is  used,  the  natural  constants  of  integra¬ 
tion  to  be  chosen  are  the  initial  peak  pressure  and  the  energy  integral. 
However,  Cole  mentions  that  it  is  difficult  to  evaluate  the  energy  inte¬ 
gral  very  accurately  from  experimental  pressure-time  curves.  Hence,  ex¬ 
perimental  data  of  peak  pressure  as  a  function  of  distance  are  often  used 

instead,  for  the  determination  of  initial  conditions  from  the  equation  for 

dP 

HR* 

I  have  not  worked  out  a  numerical  example  here  since  the  computational 
methods  are  well  known-^,  and  since  relatively  little  would  be  gained  by 
such  an  example  compared  to  the  time  required  for  computation.  (Note  the 
complexity  of  the  equations  derived  in  Appendix  XI). 

35 

OSRD  1231  (1943).  Calculation  of  detonation  pressures  for  several 
types  of  explosives. 
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COMPARISON  OF  EXPERIMENTAL  'DATA  'VITH  VALUES  CALCULATED 
■BY  THE  THEORY  OF  KIRKA'OOD  AMD  BRINKLEY 

In  Fig.  U.  I  have  plotted  the  overpressure  ratio  TT-  — — '  against 

reduced  distance  Z  r  •  The  reason  for  this  type  of  graph 

is  based  on  the  law  of  similarity^,  which  asserts  that  pressure  and 
other  properties  of  a  shock  wave  will  be  unchanged  if  the  scales  of 
length  and  time  by  which  it  is  measured  are  changed  by  the  same  factor 
as  the  dimensions  of  the  charge.  Jince  the  weight  of  the  charge  is  a 
constant  (density)  times  its  volume,  multiplying  the  weight  by  a  factor 

l/3 

f  Amounts  to  multiplying  each  linear  dimension  of  the  charge  by  f 
Hence,  the  expression  Z  for  reduced  distance  is  constructed  to  make 
all  properties  as  .functions  of  Z  independent  of  charge  weight  for  a 
given  explosive.  For  TNT,  the  density  is  such  that  it  makes  Z  -  l.Oo 
where  .  d  is  the  diameter  of  a  sphere  having  the  same  volume  as  the  charg 
and  Z  is  thus  approximately,  the  distance  from  the  explosive  source  in 

37 

charge  diameters  * 

Fig*  4*  shows  a  theoretical  curve  for  explosion  of  TNT  spheres,  as  com¬ 
puted  by  the  method  of  Kirkwood  and  Brinkley  by  the  originators  of  the 
method.  Unfortunately,  I  was  unable  to  find  data  which  are  exactly  com¬ 
parable  to  the  computed  curve. 

-*6  First  stated  by  H.  W.  Hilliar,  (British)  department  of  Scientific 
Research  and  Experiment  report  RE  142/19  ( 1 919 
37  R.  G.  Stoner  and  B.  Bleakney.  "The  Attenuation  of  spherical  shock 
waves  in  air,"  J.  Appl.  Phys.  19,.  670  -678  (July,  1948) 
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Instead,  the  following  data  are  shown  plotted  in  Fig.  4. 

Charges  Experimenters 

3  3/4  lb  pentolite  spheres  Stoner  and  Bleakney*^ 

Cylinders  of  TNT/CE,  ranging  Grime  and  Sheard^ 

between  3  and  76  lb. 

The  curve  for  the  pentolite  spheres  is  seen  to  approximate  the  theoretical 
curve  for  TNT  spheres' very  closely,  and  lies  slightly  above  it  throughout 
the  plotted  range.  This  discrepancy  may  be  explained  by  the  difference  in 
£he  explosives,  pentolite  having  approximately  a  5%  superiority  in  pressure 
level  over  TNT,  according  to  the  experimenters. 

The  plot  of  data  obtained  for  TNT/CE  is  seen  to  be  similar  in  shape  to  the 
computed  curve,  but  to  lie  a  good  deal  above  it.  This  difference  may  be 
explained  by  superiority  of  TNT/CE  to  ordinary  TNT,  and  by  the  fact  that 
the  experimental  charges  were  cylindrical  (with  diameter  equal  to  height) 
whereas  the  computed  curves  applies  to  spherical  charges.  Further,  the 
measurements  were  taken  near  the  ground,  which  may  result  in  some  distortion 
of  experimental  results. 

As  far-  as  we  may  conclude  from  the  above  discussion,  the  theory  of  Brinkley 
and  Kirkwood  thus  seems  to  be  in  good  agreement  with  experimentally  obtained 

values. 

38 
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GLOSSARY  OF  HYDRODYNAMIC  EQUATIONS 


Basic  Approaches 

Hydrodynamic  relations  may  be  expressed  in  two  different  forms,  one 
due  to  Lagrange,  the  other  due  to  Euler. 

1*  Lagrange* s  form  describes  fluid  motion  in  terms  of  the  paths  of 
individual  particles,  i.e.  it  expresses  the  coordinates  x,  y,  z, 
of  particles  as  functions  of  time  and  three  other  parameters 
(usually  taken  as  the  initial  conditions).  Lagrange's  equations 
are  usually  too  cumbersome,  but  are  advantageous  for  describing  one- 
dimensional  motion. 

2.  In  Euler's  form,  an  observer  is  considered  at  each  point  (x,  y,  z), 
and  the  happenings  at  these  points  are  observed  in  the  course  of 
time  t*  The  fluid  motion  is  described  by  giving  as  functions  of 
x,  y,  z,  t,  the  velocity  components  u,  y,  w,  of  the  particle  that 
is  at  x,  y,  z,  at  time  t. 

3*  The  relation  between  the  two  systems  is  expressed  by  the  equations 
x  x  u  (  x,  4,  a,  \) 

*  V  (x, 

i  *  w  (*;  -j,  2,-0 

where  the  dot  denotes  differentiation  with  respect  to  t  in  the 
Lagrangian  form. 

Conservation  of  Mass 


1.  Lagrange's  Form 


D  C  H. 

where  A  **  _  \  i.e.  the 

3Ui'W,c) 


Jacobian  of  x(a,  b,  c,  t),  y(a,  b,  c,  t),  z(a,  b,  c,  t);  a,b,c, 
being  the  above  mentioned  parameters. 
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2.  Euler1 s  Form 

■  .§£■*  ^ c pvj  ^ (pw)  «  ° 

This  may  be  written  more  simply  in  vector  form, 

5fL  +  PV'V/  *  O 

■at  ' 

where  the  vector  .  V/  =  ui  +  Vjj  w  ^>7  S  j  >  ^  being 
unit  vectors  in  the  x,  y,  z,  directions,  respectively. 

Conservation  of  Momentum 
1.  Lagrange  fs,  Form 

p*  *  §S‘°  ’  n  - 1?  -  f  **  If  -  ° 

which  may  be  written  a3 

*  ±-  v  p '»  o 

at  P 

where*  vr^UH+vjj+w^  *  Another  common  way  this  is 


expressed  is 


+  V  K.  *  T  VS 


3y  i 


i -  VV/1"  V/  x(  v/)  4  =  *  vs 


2#  Euler1 s  Form 


5*  =  o 


a(fu,  .  £t,u.)  .  £  ip*o  -  A *  =*  * 

s‘pv'  *  fe  Cp-v)  .  a.  <fv->  f  *  If  -  0 

^<f"'  *  W  If’*'")  *  ♦  fe  *~  5*  ■  ° 


or,  in  vector  form 


p  4* 

<  dt 

Conservation  of  Enerj 


P  -V  p(V'b)w  =  -  V 


1#  Lagrange 1 s  form 


e  +  2<j£3i  *  =  o 

r  m.  ax  ofc 

where  EV«c  =  E  *  e  bein«  internal 

energy  per  unit  mass. 
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internal  energy 
a  property,  or  1  - 
enthalpy 
constants 
functions  of  P 
over  pressure,  p-pQ 
Lagrange  coordinate 
absolute  temperature 
Shock  wave  velocity 
flow  velocity  behind  wave  front,  or  charge  weight 
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Heat  transfer  coefficient 

/* 

friction  or  viscosity  coefficient 

k 

Thickness  of 

fluijJ  disk 

TT 

~  ,  or 

Fo 

P 

density 

O 

C  *  *V 

T 

Reduced  time 

,  or  a  constant  time 

Subscript  (o)  denotes  initial  conditions  or  properties  of  the  undisturbed 
fluid. 


